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. We give a criterion for the simultaneous existence or non existence of two long-range 

Q ■ orders for two observables, at finite temperatures, for quantum lattice many body sys- 
tems. Our analysis extends previous results of G.-S. Tian limited to the ground state of 
\0 ' similar models. The proof involves an inequality of Dyson-Lieb-Simon which connects 
^ ' the Duhamel two-point function to the usual correlation function. An application to the 
OO ' special case of the Holstein model is discussed. 
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1. Introduction 

In a recent paper, G.-S. Tian derived a sufficient condition for the simultaneous ex- 
istence or non existence of long-range order of two observables in the ground state of 
quantum lattice models Here we prove a similar criterion for the case of finite tem- 
peratures. As will be seen the theorem is quite general, but for concreteness we limit 
ourselves to a specific class of fermion Hamiltonians given by (|l|) below. A concrete 
example concerning the Holstein model is given at the end and an application to the 
comparison of critical temperatures for systems with more than one order parameter is 
briefly discussed. 

Let us ffist recall the result of [|^. Let H\ be a fermion Hamiltonian of the form 

-^A ^ ^ ^xyCxa^yu ~l~ ^ ^ ^xy'^xV'y ~l" ^ ^ Jxy^x ' ^y (l) 

x,yS:A x,y£A x,y(iA 

where c^^, Cy^ are creation and annihiliation operators of spin one-half fermions on a 
(i-dimensional lattice A G Z"^, n^, and are the usual fermion density and spin at site x. 
We will always assume that the hopping matrix txy and the two body interactions Vxy, 
Jxy are integrable. This implies that -^^Tlx y£K\'^xy\ ~ ^(l) with respect to |A|, where 

(^xy ^xyi ^xyi Jxy 

Let |\E'q) be the ground state of -ffA- If is a local operator, i.e. an operator valued 
function of {c*.,?/ G V^} for a bounded neighborhood of x such that IV^I = 0(1) 
with respect to the volume |A|, we define the Fourier transform for each vector k in the 
reciprocal lattice as 

Ofc = -^ J]e-^^-^0, (2) 

V 1^1 xeA 

Suppose there exist three local observables A^, B^, such that for each a; e A 

[H^,Ax]=iiBx + uCx (3) 
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where n and u are two non-vanishing complex numbers. The resuh of [|l| states that if 
(§) and some additional assumptions (see later on) are satisfied then 

}''^.\T^^o\BtB,\^o)^0 (4) 

A^Z'i \A\ 

if and only if 

hm i-(v^^|C'+C',|M/^)^0 (5) 

A^Z'* |iV| 

Notice that if @) and (|) do not vanish they are strictly positive. In fact @ and (|) 
express the presence of long-range order for the observables Bx and C^. Indeed 

6+6, |vl.^) = E ^'"'^'"^^ <<l 0^0^ ^ E 1(^1 0^0^ lOl 

' ' ' ' x,y£A ' ' x,j/GA 

(6) 

Therefore if the left-hand side of is strictly positive one must also have 

hm hm 1(^^10+0^*^)1 >e>0 (7) 

for some e > 0. 

The result of [Q] is based on the following inequality which follows from the commutator 
relation (^). For any k 

■{^^\[Kt,[H,Kk]]\^^y/' (8) 

where Kk = [HA,Ak]- If the right hand side of (^) is 0(1) uniformly in |A|, the two 
terms in the square are necessarily of the same order. More precisely assuming that the 
right-hand side is 0(1), if one of the correlation functions is 0(|A|), this is also the case 
for the other. The reader is referred to for several concrete examples. Let us note 
that if u vanishes (resp. fi vanishes) (P) implies that (^I/q | B^Bk |\I'o ) = 0(1) for all k 
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which means there is no long-range order for (resp. C^) 0. In the next section we 
prove a similar criterion for the case of finite temperatures. 

2. Finite temperatures 

Let {Ox) A = {OxQ~^^^^~^^'^^ ) / Za where the trace is over the usual Fock space, 
A^A = Tlx&K -^A ~ '^^ e~f^^ is the partition function at inverse temperatures (3. 

The generalisation of (§) to finite temperature involves the Duhamel two-point function 
whose definition we recall here. The Duhamel two-point function of two operators F 
and G is defined by 

Za Jo 

It is symmetric (F, G)a = {G, F)a and satisfies a Cauchy-Schwarz inequahty 

|(F,G')a|2<(F+,F)a(G+,G)a (10) 
We also recall the following useful identity 

{[F,G])A = i[F,(3HA],G)A (11) 
The following lemma generalizes the inequality (|[). 

Lemma 

Let Ax, Bx and Cx be three local observables satisfying (Qj then for all k 

{\li\^{Bt,B,)A-\u\\l{Ct.C,)^ <-^{[At,k,])A (12) 

Proof: 

Setting F = i+ and G = [Ha, ifc] = Kk in (|Tl|) 

h[Al, [Ha,Au]])a = {[AI^HaI [Ha,Aj,])a = {[Ha^Au^, [Ha,Au])a (13) 
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Replacing the Fourier transform of (|^) in the right-hand side of ([T^) we obtain 



\l2\\B^, Bk)A + W\\C^, Ck)A + i^*iy{B^, Ck)A + iy*KC^, 4)a (14) 



We can estimate the last two terms on the right-hand side of (|Tj) using the Cauchy- 
Schwarz inequality 

\fi*u{B+,Ck)A + 'y*fi{Cl,Bk)A\ < \H{\{B^,Ck)A\ + \iC^,Bk)A\) 



<2\^^v\^{Bt.B,)A\l{Ct.C,)A (15) 
Noticing that the term /i*z/(_B^, CtjA + ^*fJ'{C^, Bk)k is in fact real, we have 



^l*v{Bt.Ck)^ + y*^l{Ct.Bk)A > -2\^^v\^J {B^ .B,,)a\I {Ct .Ck)A (16) 
From ([T^) and (|16]) we get 

> {\^^\^{Bt,Bu)A-\u\^{Ct.C,)^ (17) 

which is the desired inequality. To extract information about the correlation functions 
of B^ and we have to connect these to the Duhamel two-point function. This will be 
done through upper and lower bounds on the Duhamel two-point function which involves 
only the symmetrized correlation function because the Duhamel two-point function is 
symmetric. The upper bound follows from convexity 

{F+,F)a<\{F+F + FF+)a (18) 

The lower bound was found by Dyson, Lieb and Simon 

(F+, F)a > \{F^F + FF+)Af{hA{F)) (19) 
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where 



and the function f{u) is defined imphcitely for u > by 

/(utanhn) = — tanhu (21) 
u 

The function f{u) is continuous convex and strictly decreasing with hm„^o f{u) = l and 
hm^,_oo f{u) = 0. 

For any local observable Ox, we set 

OA{k) = ^-{dtd, + d,dt)A (22) 



From (]12|) we get the two inequalities 



|/.|(4+, 5,)f > W\iCt,C,)T - -^{[At,K,])Y' (23) 
W\iC^,C,)T > MB^,B,)]l" - ^{[^t,Kk])T (24) 



Using the upper and lower bounds ([TSD and (0) we have 

If^mk)'/' > WlCAiky/'fihAiC,))'/' - -^{[At,k,])Y' (25) 
\u\C^iky/' > \^\B^{kf'^f\h^{B,)Yl' - (26) 

From the two last inequalities and the behavior of f{u) for m — > 0, it is easy to deduce 
the following 

Theorem 

Assume there exist three local observables A^, B^, Cx satisfying (Qj. Suppose also that 
for a given k, there exist three positive constants a^, bk and Ck independent of |A| such 
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that [H^,Ak]])A < ak, [H^,Bk]])A < bk and {[C+ , [HA,Ck]])A < Ck. Then for 

any temperature f3^^ 

Ita^B.W^O (27) 



if and only if 



lim -^CJk)^0 (28) 



In the last theorem, it is important that the expectation values of the double com- 
mutators are 0(1) with respect to the volume. In the Appendix, we show that this is 
indeed the case for the Hamiltonian (|l]) and a large class of observables. It turns out 
however that this assumption can often be checked explicitely in concrete applications. 

3. Applications 

We wish to present an illustration of the theorem for the dynamic Holstein model. This 
is a model for molecular crystals which retains the interaction of itinerant electrons with 
breathing modes of the molecules. These modes are described by Einstein oscillators of 
frequency u attached to each site x G A. Let Qx, Px be the position and momentum of 
the oscillators (with [qx,Py] = iSxy)- The Hamiltonian of the model is 

Choosing = Px we obtain 

[Ha, Px] = ^U{nx + iuJ^Qx (30) 

so that Bx = {ux — \)-, Cx = Qxi fJ' = iU and u = iu^. To apply the theorem we have to 
check that {[O^, [H\,Ok]])A = 0{1) for Ok = A^, Bk and Ck- By exphcit calculation we 
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have 

[i+,[ifA,ifc]]=^2 (31) 

[4+, [H^, B,]] =i- J2 2(cos(A; ■ {x - y)) - l%yctcy (32) 

[C'+,[ifA,d]]=l (33) 

The double commutators with and Ck are constant and thus equal to their expectation 
values. For Bj. we have 

[ifA,5fe]])i < ^ E l2(cos(^ ■ - y)) - 1)1 1^^.11 (4c.) Ia < m E I^-j'I = 

' ' x,yeA ' ' x,yeA 

(34) 

where we used that for fermionic operators \{c'^Cy)\ < 1. Thus we can conclude that 
limA_^2d jx\{in^k - l)ink - |))a 7^ if and only if limA-^z^ |i^(g-fc^fc)A 0. 

It follows from this result that if a Peierls instability occurs at a given k (i.e. a 
distortion of the lattice with long-range order for q^) there is at the same time long- 
range order in the electronic density. The occurence of a Peierls instability for k = ti 
at low temperature and half-filling has in fact been proven for the static model where 
the kinetic energy of the phonons in (^) is absent [0. The present analysis remains 
unchanged in the static case because even if the term p^. is absent from the Hamiltonian 
we still have the right to choose = Px so that (|30D remains unchanged. 

In one dimension, for the spinless case the decay of the fermionic correlation functions 
for the ground state has been analyzed rigorously [^] at weak coupling U by using 
renormalisation group techniques. The results of |^ imply that these correlations decay 
algebraically and therefore we can conclude that no long-range order is present at the 
level of as well. 

An interesting application of our theorem concerns the comparison of critical temper- 
atures in systems with several types of degrees of freedom and several order parameters. 
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An example is the Holstein model itself. Another example are the so-called ferrofluids 
where one has magnetic (spin) degrees of freedom as well as translational degrees of 
freedom (see for example P). A relevant question is whether there exist two different 
critical temperatures associated to each order parameters or only one critical tempera- 
ture where both order parameters simultaneously develop a non zero value. For example 
in the ferrofluids, it is of interest to decide whether liquid-vapor and ferro-paramagnetic 
phase transitions occur at the same critical temperature or not. Although our theorem is 
a statement about existence of long-range order rather than about the order parameter, 
it suggests that when it is applicable there is only one critical temperature. This fact 
had been proven earlier for the particular case of the static Holstein model in 0. 

A related question that our basic result leaves open is whether the behavior of the two 
correlation functions at a critical point is governed by two independent length scales. 
Our discussion above suggests that (when our result is applicable) this is not the case 
and there is a single diverging length scale involved in both correlations. 
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In order to apply our theorem, it is necessary that all the expectation values of the 
double commutators are 0(1). In this appendix, we check that this is the case for the 
Hamiltonian (|l|) and a large class of observables. The double commutators we have to 
discuss are 



Appendix 




(35) 
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where is a local observable centered around z of the form 
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Oz — • • • C^+5j^C2+e-^ . . . Cz+eM (36) 

where i = 1, . . . ,N and e^, j = 1, . . . , M are fixed vectors and N and M are fixed 
numbers. The subsequent discussion provides a proof that (35) is 0(1) when + M is 
even. If + M is odd we have only heuristic arguments presented at the end. 

We give the detailed calculation for the kinetic term of (|l]) (for simplicity we omit 
spin indices). The interaction terms can be treated in a similar way. Let us evaluate the 
commutator of Oz with the kinetic term 

N 

[ txyC^Cy, Oz] = t^ySy^z+SiC^+Si ' ' ' '^2+<5,_ 1 '^"T '^z+5j+ 1 • • • ^t+Sff^z+ei ■ ■ 'Cz+eM 

x,y(iK i=l x,y£A 

M 

— txySx^z+ejC'^+g_^ . . . C'^^g^Cz+ei ■ ■ ■ Cz+tj^xCyCz+ej+-i ■ ■ ■ Cz+€m (^7) 

j=l x,y£A 

The commutator of O^, with the first term of the right-hand side of the last equation 
reads 



N 

i=l x,y£A 



C'Cl 



h5i+i • • • '^z+S^'^z+ei ■ ■ ■ Cz+e 



^z+5i ■ ■ ■ ^z+5i-i^x ^z- 



-S,+i ■ ■ ■ Ct+Sj^'^z+'^I ■ ■ ■ Cz+eM<^t'+eM ■ ■ ■ '^t'+ei'^z'+SN ■ ■ ■ <^z'+5i) (38) 



Using the fermionic commutation relations, we can drive in both terms the creation 
operators to the left and the annihiliation operators to the right. Doing this we obtain 
a finite number of terms involving two Kronecker symbols of the type Sy^z+5,Sz'+ii,z or 
Sy,z+SiSz'+ti,x, with /i a linear combination of 5's and e's, plus a supplementary term 
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containing only one Kronecker symbol. The latter one is 



N 

(-i)^^i - J2 E t^y^,^+^M'+.M ■ ■ ■ ■ ■ ■ 



1=1 x,yeA 

ct+S,+ i ■ ■ ■ ct+Sj^Cz'+Si^ ■ ■ ■ Cz'+5iCz+ti ■ ■ ■ Cz+eM) (39) 



We can now replace these expressions into (^5|) . Let us first suppose that M + is 
even. In this case the term (|39D vanishes and only the terms with two Kronecker symbols 
contribute. It follows from the previous discussion that the contribution of (|38| ) to (^) 
is bounded by 

1 ^ 

mE E M^,y,z,z')\t,y\\{...c^...c...)A\ (40) 

' ' i=l x,y,z,z'eA 

where Aj(x, y, z, z') is one of the following: 5y^z+5i5z'+^l,z or Sy^z+Si^z'+^i^x- The expectation 
value of the product of c* can be bounded by 1 since the norm of the fermionic operators 
is less than 1. Then using the integrability of t^y, the two Kronecker symbols and the 
factor 1/1 A|, we conclude that (^) is bounded by some constant uniformly in the volume. 
It is possible to perform a similar discussion for the second term in the right-hand side 



of (0) to conclude that the double commutators of local observables of the type (^) 
with H have to be 0(1) with respect to the volume in the case N + M even. 

In the case where + M is odd, the discussion for the terms with two Kroneckers 
symbols is the same as in the case of iV + M even. However, there remains the term 
(^). Replacing it in (|35|) produces a term bounded by 



N 



' ' i=l x,y,z,z'&A 

Cz'+Sn ■ ■ ■ Cz'+5iCz+ei ■ ■ ■ Cz+eM) a\ (^1) 

In this last expression, the averages of the product of creation and annihiliation operators 
are matrix elements of the reduced density matrix Pm+n- In Yang conjectured that 
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long range correlations (ODLRO) may appear in such matrix elements only for N + M 
even. We are not aware of a general proof of this conjecture but it suggests that the 
matrix element in (41) decays as |x — z'| ^ oo. Assuming the decay is quick enough the 
sum over z' in (41) converges, which then leads to an upper bound 

iV 

r^Xl 5Z ^y,^+s^\txy\ (42) 

1=1 x,y,z£A 

where c is a constant independent of the volume. Summing over z, using the integrability 
of the hopping amplitudes and the term 1/|A|, we have that (^2]), and thus (0), are 
bounded by a constant uniformly in the volume. 
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